Quantum freeze of fidelity decay for chaotic dynamics 
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We show that the mechanism of quantum freeze of fidelity decay for perturbations with zero 
time-average, recently discovered for a specific case of integrable dynamics [New J. Phys. 5 (2003) 
109], can be generalized to arbitrary quantum dynamics. We work out explicitly the case of chaotic 
classical counterpart, for which we find semi-classical expressions for the value and the range of the 
plateau of fidelity. After the plateau ends, we find explicit expressions for the asymptotic decay, 
which can be exponential or Gaussian depending on the ratio of the Heisenberg time to the decay 
time. Arbitrary initial states can be considered, e.g. we discuss coherent states and random states. 
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The question of stability of quantum time evolution 
with respect to small changes in the Hamiltonian has re- 
cently attracted lot of attention 0, Q • This question is 
particularly important in the context of quantum infor- 
mation processing The central quantity for describ- 
ing quantum stability is the fidelity F(t) — | (-(/>(£) ^(i))! 2 
where \ip(t)) = U {t)\vb) and |^(*)> = U s (t)\ip) are un- 
perturbed and perturbed time evolutions, of perturba- 
tion strength S, respectively, starting from the same ini- 
tial state \ip). Let the evolution operator be written as 
time-ordered product U$(t) = Texp(— i j Q dt' H$(t')/K) 
in terms of (generally time-dependent) Hamiltonian 
Hg(t) = Ho(t) + 6H'(t). In this Letter we either as- 
sume that Hs(t) is autonomous (time-independent), or 
more generally, periodically time dependent with some 
period r, Hg(t + t) = H$(t). Then the time is mea- 
sured in discrete units of r, namely t = tit, and the 
former (autonomous) case is simply obtained as the limit 
t — ► 0. The perturbed propagator for one time step can 
be written as U§(t) = Uq(t) exp(— iVrS/H) in terms of 
a hermitean perturbation V which in the leading order 
perturbs the Hamiltonian V = H ' + 0(tS). 

It has been shown Q that for classically chaotic sys- 
tems and for sufficiently strong perturbation and coher- 
ent initial state the fidelity decay is given by classical 
Lyapunov exponents, and this phenomenon has been re- 
cently explained solely on the basis of classical dynamics 
0. On the other hand, for sufficiently small 5, one can 
express fidelity decay in terms of a power series in S where 
coefficients are given as time-correlation function of the 
perturbation [I| . Using this approach one can derive uni- 
versal forms of fidelity decay in both cases of classically 
regular and chaotic dynamics and express all time-scales 
solely in terms of classical quantities and h. 

The starting point of our analysis is the representation 
of fidelity F^ n ' — F(tvt) in terms of expectation value |5j 



F (n) = | / (n)|2 j /(") = (^|M 5 ( "V) (1) 



of the echo operator Mg := Uo{—nT)Us{nT) which is 



the propagator in the interaction picture. Namely 



(n) 



Texp ( -i^E„(V) ) , 



(2) 



where, for any operator A, T, n (A) :— TY2n'=o^ n ' an< ^ 
An := Uq(— nT)AUo(nT). In case of continuous time: 
M s (t) = f exp(-i££(V,i)), with T,{A,t) := f*dt'A(t'), 
A(t) := U (-t)AU (t). Approach @ using the power law 
expansion of in 5 gives to the second order 

F {n) = 1 - ^ {(^(F)) - (S„(P)) 2 } + 0{5 A ) (3) 

where (•) := (tp\ • \vj). Equivalently, this useful formula 
can be expressed in terms of time-correlation function 
C(n>") = {V n ,V n „) - (V n .) (V n »), namely = 1 - 

(ts/h) 2 E^'~n"=o n ") + ■ ■ ■ The rule of thumb sa y s 

that slower decay of correlations, i.e. stronger fluctuation 
of £„(V), imply faster decay of fidelity, and vice versa. 

Particularly interesting special situation arises when a 
time averaged perturbation V := lim„^ 00 (rir) _1 S„(U) 
equals zero. In general, the perturbation can be decom- 
posed into the diagonal and residual part V = V + V TCS . 
The part V which commutes with the unperturbed evo- 
lution Uq and is thus diagonal in its eigenbasis can some- 
times be put together with the unperturbed Hamiltonian 
Hq. This is customary in various quantum mean field ap- 
proaches. It is thus interesting to question the stability of 
quantum dynamics with respect to residual perturbation 
only (i.e. when its diagonal part exactly vanishes V = 0). 
This problem has been addressed for the particular case 
of perturbed integrable dynamics and very interesting 
results on extreme stability of quantum dynamics have 
been found termed as 'quantum freeze of fidelity'. 

In this Letter we show that the phenomenon of quan- 
tum freeze, namely the saturation of fidelity to a plateau 
of high value ^is much more general and robust as it ap- 
pears in Ref.jy], and applies to arbitrary quantum evo- 
lution provided only that the perturbation is residual, 
V = 0. In particular, we work out in detail the important 
case of dynamics with fully chaotic classical counterpart. 
We compute the plateau value (scaling as 1 — const(5 2 
within the second order), its range scaling as 1/5, and 
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the rate of the asymptotic decay after the plateau ends 
(which is cither gaussian or exponential), quantitatively 
in terms of the underlying classical dynamics, effective 
perturbation 6 and the effective value of Planck constant 
h. The phenomenon may find useful application in quan- 
tum computation where fidelity error is predicted to be 
very small and frozen in time (for sufficiently small <5) 
provided only that the diagonal part of the error in each 
gate can be cured by some other means. 

In the autonomous case (r — ► 0), provided that the 
spectrum of H is non-degenerate (which is true for a 
generic non-integrable system), the perturbation is resid- 
ual iff it can be written as a time derivative of some ob- 
servable W i.e. a commutator with H , V — i[Ho, W] — 
(d/dt)W. Generalizing to the discrete, time-periodic case 
we shall assume that the perturbation is of the form 

V = - (Wi - W ) = - (Uo(-t)WUo(t) - W) . (4) 

T T 



We shall now apply the Baker-Campbell-Hausdorff ex- 
pansion e A e B = exp(A + B + (1/2) [A, B] + . . .) to the 
time-ordered product J2J) and rewrite the echo-operator 



M] n) = exp <J ~ ( E n (V)5 + ^T n S 



1 



(5) 



where T n := if YZ>=o E™»i„'[K/, V n »]. It is interest- 
ing to note 6] that all matrix elements of r„ grow with 
n not faster than oc n [jj . This becomes obvious for the 
special form of perturbation Q for which it follows 



£„00 = W n -Wo 

Tn 



(6) 



E n (iZ)-i[W ,Wn], i2:=-^[W ,Wi], (7) 



so the operator T n is also a time sum/integral of a time- 
dependent operator R, minus a sort of time-correlation 
function which shall be neglected for systems with strong 
decay of correlations studied below. In the continuous 
time case, R = j-[W, (d/ dt)W] = h- 2 [W, [W, H ]} and 
r(t) = J*dt'R(t') - j-[W(0),W{t)]. We note that, pro- 
vided W has a well defined classical limit H — > 0, then 
also V, R and r„ have well defined limits since ¥[•>•] 
can be replaced by a Poisson bracket. This is what we 
shall assume below, as well as that the limiting classical 
dynamics of Uq, is fully chaotic. 

Comparing the two terms in the BCH exponential © 
we note that there should exist a time-scale t 2 ~ <5 _1 , 
such that if nr < t 2 then the first term T, n S dominates 
the second one ^T n S 2 (and higher 0). So, let us first 
consider the case nr < t 2 . Then we can neglect the 
second term and write the fidelity amplitude 



/W = (exp(-i(W n -W )S/h)) 



(8) 



Expanding f( n > to the second order in S, we find = 

1 - g («§ + < - C n - C*) where k\ := (W 2 k ) - (W k )\ 
C n := (W n W Q ) - (W n ) (Wo). Using Cauchy-Schwartz 



inequality |C„| < kqk u and the fact that for a bounded 
operator W the sequence K n is bounded, say by r, we find 
a freeze of fidelity 1 — < 4j-?r 2 , nr < t 2 oc S^ 1 , for 
arbitrary quantum dynamics, irrespective of the existence 
and the nature of the classical limit. 

Let us further assume that, due to mixing property 
of classically chaotic dynamics, time correlations van- 
ish semiclassically beyond some mixing time-scale t\, 
C n — * 0(h) for nr > t\, and quantum expectation values 
become time-independent and equal, to leading order in 
h, to the classical averages over an appropriate invari- 
ant set (A) cl :— f dfiA c \. Hence, between t\ and ti, the 
fidelity freezes to a constant value @ 

Fplat + 4) , «d : = (^ 2 )cl - (^)cl • (9) 

Considering two interesting extreme examples of initial 
states, namely coherent initial states (CIS) and random 
initial states (RIS) we find: For CIS Kq oc ft can be ne- 
glected with respect to n 2 h whereas for RIS K n does not 
depend on time hence ^cl" 

So within the linear 
response approximation 1 — -Fpi a t is universally twice as 
large for RIS than for CIS. It is also worth to stress that 
quantum relaxation time t\ oc log H for CIS while t\ oc h° 
is simply the classical mixing time for RIS. 

One can go beyond the linear response in approxi- 
mating (jSJ) using a simple fact that in the leading or- 
der in h quantum observables commute, and as before, 
that for nr > t\ the time correlations vanish, namely 
(exp(-if (W„ - W ))) w (exp(-ifwy) (exp(i$W )>: 



Fput » I (exp (-iWS/H)) cl (exp {iW S/h))\' 



(10) 



Defining a generating function in terms of the classical 
observable W c i, G(z) :— (exp(— izW)) cl , one can com- 
pactly write w \G(S/h)\ 2 for CIS (neglecting local- 
ized initial state average with Wo) and ~ \G(S/h)\ 4 
for RIS, satisfying universal relation F^ » (i^^t) 2 - 
Curiously, the same relation is satisfied for the case of 
regular dynamics 0. If the argument z = S/h is large, 
the analytic function G(z) can be calculated generally by 
the method of stationary phase. In the simplest case of 
a single isolated stationary point x* in N dimensions: 



\G{z)\ 



N/2 



detd Xj d Xk W(x*)\ 



-1/2 



(11) 



This expression gives an asymptotic power law decay of 
the plateau height independent of the perturbation de- 
tails. Note that for a finite phase space we will have 
oscillatory diffraction corrections to eq. (|llf> due to a fi- 
nite range of integration J d/i which in turn causes an 
interesting situation for specific values of z, namely that 
by increasing the perturbation strength 8 we can actually 
increase the value of the plateau. 

Next we shall consider the regime of long times nr > 
ti. Then the second term in the exponential of JSJ dom- 
inates the first one, however even the first term may not 
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be negligible. Up to terms of order 0(nS 3 ) we can factor- 



ize eq. © as M { 5 n) « exp(-i|(W„ - W )) exp(-i£r„). 

When computing the expectation value = (M$ ) 
we again use the fact that in the leading semiclassical 
order the operator ordering is irrelevant and that, since 
tit 3> t\, any time-correlation can be factorized, so also 
the second term of r„ (0 vanishes. Thus we have 



« F, 



plat 



exp 



S 2 , ^ 

-i Th Z n (R) 



tit > t 2 . (12) 



This result is quite intriguing. It tells us that apart 
from a pre-factor .Fpiat, the decay of fidelity with resid- 
ual perturbation is formally the same as fidelity decay 
with a generic non-residual perturbation, eqs. IIIL'I) . 
when one substitutes the operator V with R and the 
perturbation strength 8 with 8r = 5 2 /2. The fact 
that time-ordering is absent in eq. i|12[) as compared 
with @ is semiclassically irrelevant. Thus we can di- 
rectly apply the general semiclassical theory of fidelity 
decay [|| , using a renormalized perturbation R of renor- 
malized strength Sr. Here we simply rewrite the key 
results in the 'non-Lyapunov' perturbation-dependent 
regime, 5rt < h. Using a classical transport rate 



( r:=lim I ^ 0?5 i 7 «S2( J R)) cl 
an exponential decay 



(E„(iZ)) cl ) we have either 



1 r— 
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FIG. 1: F(t) for the kicked top, with 8 = 10" 3 (a), and 
8 = 1CP 2 (b). In each plot the upper curve is for CIS and the 
lower for RIS. Horizontal chain lines are theoretical plateau 
values (1101 . vertical chain lines are theoretical values of ti 
l|15|l . The full circles represent calculation of the correspond- 
ing classical fidelity for CIS which follows quantum fidelity up 
to the Ehrenfest (log K) barrier and exhibits no freezing. 



F {n> « F plat exp ( -7^ anT ) > nT < *H 



or a (perturbative) Gaussian decay 
plat6P l 2h 2 t H 



tit > in- 



(13) 



(14) 



tn = tN '/(2s) is the Heisenberg time, where N ~ ftr d 
(in d degrees of freedom) is the total dimension of the 
Hilbert space supporting the time evolution and s is the 
number of different symmetry classes (of possible discrete 
symmetries of H ) carrying the initial state \ip). This is 
just the time when the integrated correlation function of 
R n becomes dominated by quantum fluctuation. Com- 
paring the two factors in I|13I14|) . i.e. the fluctuations of 
two terms in (JSJ , we obtain a semiclassical estimate of t 2 




8 ' aS 2 



(15) 



Interestingly, the exponential regime (|13f) can only take 
place if ti < £h- If one wants to keep -Fpiat ~ 1, or have 
exponential decay in the full range until F ~ 1/7V, this 
implies a condition on dimensionality: d > 2. The quan- 
tum fidelity and its plateau values have been expressed 
(in the leading order in h) in terms of classical quantities 
only. While the prefactor -F p i a t depends on the details 
of initial state, the exp- factors of I|13I14|) do not. Yet, 
the freezing of fidelity is a purely quantum phenomenon. 



The corresponding classical fidelity (defined in |5j) does 
not exhibit freezing. Let us now demonstrate our theory 
by numerical examples. 

First we consider a quantized kicked top as an ex- 
ample of one-dimensional system (d = 1). The sys- 
tem is described by quantum angular momentum J Xj y iZ 
with (half)integer modulus J and the one-step propaga- 
tor U = exp (—ia>J 2 /2J) exp (— mJ y /2). We have chosen 
a = 30 ensuring fully chaotic corresponding classical dy- 
namics, with angular momentum coordinates on a unit 
sphere named as x,y,z, and J = 1000 determining the 
effective Planck constant h = 1/J = 10~ 3 . The per- 
turbation is chosen as V — (J 2 — J 2 )/2J 2 associated 
with W = J 2 /2J 2 . The initial state is either RIS (with 
Gaussian random expansion coefficients) or SU(2) CIS 
centered at (ip, 6) = (1,1). In both cases the initial 
state is projected on an invariant subspace of dimension 
N = J spanned by Hqe = {|2m) - | -2m), \2m- 1) + 
| — (2m— 1)); m — 1, . . . , J/2} where \m) is an eigenstate 
of J z . We first checked the plateau. Within the linear 
response © we have to evaluate only = 1/45 for 
the corresponding classical observable W c \ = z 2 /2, giv- 



1 - {8J) 2 k% F™ = 1 - (5J) 2 2n 2 cl . These 



values give good agreement with the fidelity for weak per- 
turbation 8 = 10 -3 shown in fig. whereas for strong 
perturbation 8 = 10~ 2 shown in fig. \I]p the theoretical 
values ljTU)l of F p i at , expressed in terms of the generating 
function G(z) for CIS/RIS, have to be calculated exactly, 
and indeed the agreement is excellent. Integration over 
the sphere yields G(SJ) = v^erf (e i7r / 4 ^/SJ/2). Com- 
paring with the asymptotic general formula for G(z) l|ll|) 
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FIG. 2: Long-time Gaussian decay for CIS of a single kicked 
top for the same parameters as in fig.^,. Full curve is a direct 
numerical evaluation, empty circles are numerical calculation 
using the renormalized strength 8r and operator 7?, while the 
chain curve gives the theoretical decay 1141 . 
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FIG. 3: Long-time fidelity decay in two coupled kicked tops. 
For strong perturbation 8 = 7.5 ■ 10~ 2 (a) we obtain an ex- 
ponential decay, and for smaller 8 = 2- 10~ 2 (b) we have a 
Gaussian decay. Meaning of the curves is the same as in fig.|2l 

we now also find a diffractive contribution due to oscilla- 
tory behavior of the complex erf-function. Small (quan- 
tum) fluctuations around the theoretical plateau values 
in fig.^lie beyond the leading order semiclassical descrip- 
tion. In fig. ^we also demonstrate that the semiclassical 
formula (|15fl for i 2 works very well. Long-time Gaussian 



decay for the parameters of fig. is shown in fig. [2] 
Here we compare a direct numerical calculation with the 
numerical calculation using a renormalized perturbation 
strength 5r and the effective perturbation operator (JJJ 
R = jjj(J x JyJz + JzJyJx), and with the theoretical pre- 
diction l|14|) where the classical dynamics of R c \ = —xyz 
gives a = 5.1 • 10 -3 . 

To demonstrate the possibility of clean exponential 
long-time decay of fidelity \Y6\ we look at a system of 
two {d = 2) coupled tops J\ and J 2 given by a propaga- 
tor U = exp (— ieJ z i J z2 ) exp (— i7rJ y i/2) exp (— i7rJ y 2/2), 
with the perturbation generated by W = ® 1 + 1 ® A 2 , 
where A = J^/2J 2 for each top. We set J = l/h = 100, 
and £ = 20 in order to be in a fully chaotic regime. 
The initial state is always a direct product of SU(2) co- 
herent states centered at ((p,0) = (1,1) which is sub- 
sequently projected on an invariant subspace of dimen- 
sion TV = J(J + 1) spanned by {Hoe ® W r } sym , where 
H r = H\ 7Yoe and {-}sym is a subspace symmetric with 
respect to the exchange of the two tops. The results of 
numerical simulation are shown in fig. [3J Here we show 
only a long-time decay, as the situation in the plateau 
is qualitatively the same as for d — 1. For large enough 
perturbation one obtains an exponential decay shown in 
fig. 13 (a), while for smaller perturbation we have a Gaus- 
sian decay shown in fig. [3] (b) . Numerical data have been 
successfully compared with the theory i|13ll4fl using clas- 
sically calculated a — 9.2-10 -3 , and with the "renormal- 
ized" numerics using the operator R Q. 

In this Letter we discussed a freeze of fidelity for arbi- 
trary quantum evolution provided only that the diagonal 
part of the perturbation in the basis of the unperturbed 
evolution exactly vanishes. The value of the plateau can 
be arbitrary close to 1 and can span over arbitrary long 
time-ranges for sufficiently small strength of perturba- 
tion. We worked out in detail the case of systems with 
fully chaotic classical limit. Our result is predicted to 
have immediate application to quantum information pro- 
cessing. If combined with the inequality between the pu- 
rity I(t) of a reduced density matrix of a bipartite quan- 
tum system and the fidelity, namely J(t) > |-F(i)| 2 we 
predict that decoherence as characterized by I(t) should 
also exhibit freeze for the particular class of perturba- 
tions. Useful discussions with T.H.Seligman and finan- 
cial support by the Ministry of Education, Science and 
Sports of Slovenia and DAAD19-02-1-0086, ARO United 
States are gratefully acknowledged. 
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